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Abstract 

Analytical expressions for the distribution of the length of chords 
corresponding to the affine invariant measure on the set of chords are 
given for convex polygons. These analytical expressions are a computa¬ 
tional improvement over other expressions published in 2009 and 2011. 

The correlation function of convex polygons can be computed from the 
results obtained in this work, because it is determined by the distribu¬ 
tion of chords. 

An analytical expression for the multi-chord distribution of the length 
of chords corresponding to the affine invariant measure on the set of 
chords is found for non convex polygons. In addition we give an al¬ 
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pression. The results also apply to non simply connected polygons. 

Keywords: Chord length distribution function; polygons; diagonal pentagon; 

correlation function 

2010 Mathematics Subject Classification: Primary 60D05; 52A22 


*E-mail: r.garcia-pelayo@upm.es 


1 



1 Introduction 


The distribution of the length of chords of a plane figure is a topic that became 
popular in connection with Bertrand’s paradox in 1907 [30j [12]. The study of 
this paradox led mathematicians to realize that different problems lead to 
different measures on the set of chords. The relations between these measures 
on the set of chords of convex sets of the Euclidean space have been studied 
[3D El El]. Among these measures there is only one (up to a factor), usually 
denoted by //, which is invariant under translations and rotations (see the 
first three chapters of [T5] or chapter 1 of [30] or [38]). The distribution of 
the length of chords corresponding to /r (and other measures) had been found 
analytically for the equilateral triangle [H], the general triangle [I6l III] , the 
circle, the rectangle [13], the ellipse [35], the regular pentagon [2], the regular 
hexagon (22]. the trapeze rectangle ra and the isosceles trapezium [39], to 
name only plane figures. In the last years the distribution for regular polygons 
has been found [281 5]. In 2009 Ciccaricllo gave an analytical expression for 
the correlation function of convex polygons [10] and showed that its second 
derivative is the chord length distribution (Appendix C of [10]. In 2011 the 
correlation function of convex polygons was found by H. S. Harutyunyan and 
V. K. Ohanyan [29] differently. First they found the distribution for each pair 
of sides of a convex polygon following an approach based on Pleijel’s identity 
( 0 , page 156). Then they wrote the distribution for convex polygons as the 
sum of the distributions for each pair of sides. The work of these three authors, 
however, doesn’t furnish an efficient tool to compute distributions of the length 
of chords. For example in section III of [ID] and in section 5 of (29] one can 
see that it takes them a lot of work to compute the distribution of the chords’ 
length for an isosceles right triangle and a rhombus, respectively. 

In section 2 of this article we find the distribution for triangles by direct 
geometrical arguments. The equilateral case was found in 1961 w, page 
57) and the general case was found in 2010 by Ciccariello [IT]. However, 
his final expressions are much more involved than the ones presented here. 
Ciccaricllo has checked on an acute and an obtuse triangle that both methods 
match [9]. In section 3 we find the distribution of the length of chords of two 
concurrent segments, as a step towards finding the distribution for a pair of 
sides in section 4. That the latter distribution can be derived from the former 
has been acknowledged earlier [6]. In order to find the distribution of the 
length of chords for two concurrent segments, we have used the distribution 
for triangles as a basis, instead of Ciccariello’s approach [10] or Pleijel’s identity 
( 0 , page 156). This and the use of the indicator function of an interval yields 
mathematical expressions which are computationally more friendly then the 
corresponding ones in [29]. In section 5 we also write the distribution for 
convex polygons as the sum of the distributions for each pair of sides and, in 
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addition, we prepare the analytical expressions for their input to be Cartesian 
coordinates. The graph of the distribution of the length of chords for a rhombus 
can be produced in less than a second on a personal computer by simply giving 
the coordinates of its vertices. In section 6 we compute the distribution of the 
chords’ length of the Mallows and Clark dodecagon [33], which has theoretical 
and practical interest (see the Introduction of |24j). 

While much is known about the distribution of random chords in convex 
bodies [36], the distribution of the length of chords in non convex figures has 
features which are absent in the convex case. In fact, one can tell whether a 
body is convex or not from its distribution of chords (HU. p. 61, proposition 
4). This is also apparent when one tries to find the distribution by Monte 
Carlo methods [33]. Despite these difficulties the distribution of the length of 
chords of several non convex figures has been computed or found analytically 
[221 [2S1 [D a ESI uni EZl S0U2S]- A distribution for the length of chords can be 
found by taking the second derivative of the correlation function HB1S21HD]- 
This distribution, however, is a signed distribution for non convex bodies, 
that is, the probability may take negative values. This happens for chords for 
which one end lies inside the body and the other end lies outside the body. 
Considerations of such chords follows naturally from the interpretation of the 
second derivative of the correlation function. When we restrict ourselves to 
chords which lie fully within the body one may, for concave bodies, consider 
the sum of the lengths of the interceptions of a straight line by the body as 
the random variable, and this is the One-Chord Distribution (OCD) [21], 

Here we take what from a human point of view (perhaps not from a math¬ 
ematical point of view) is the most natural choice: we consider only chords 
which lie fully within the body and consider each segment independently from 
others which may lie on the same straight line. This is called the Multi-Chord 
Distribution (MCD) [21] . Here we find an analytical expression for the MCD 
corresponding to the invariant measure /i for non convex polygons. We also 
give an algorithm to find the MCD which is numerically more efficient than the 
said analytical expression and use it to compute the MCD for the five-pointed 
star. The results also apply to non simply connected polygons. 

In sections 2-6 we treat the case of convex polygons and in sections 7-14 
the case of concave polygons. We provide a link where two Mathematica hies 
with implementations of the procedures used in the convex and concave cases, 
respectively, can be downloaded. [ 19] . 

The distribution of chords’ length has applications in the analysis of images 
provided by a microscope [7] and others (see, for example, the first paragraph 
of |2Tj). The importance of the correlation function of polygons in small angle 
scattering has been the subject of recent research m hd]. The correlation 
function of convex polygons can be computed from the results obtained in this 
work, because the distribution of chords determines the correlation function 
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m, for a clear exposition of this in English see Appendix A of [32]). 


2 Distribution of the length of chords in a tri¬ 
angle 

In this work we shall use repeatedly these two facts: the affine invariant mea¬ 
sure p of a set of parallel straight lines is proportional to its cross section and 
the affine invariant measure p of a sheaf of straight lines is proportional to the 
angle that they span. See, for example, references |30, 32] [38] for details. We 
prefer to find the distribution of measure rather than the distribution of prob¬ 
ability to avoid dealing with normalization factors. It is better to normalize, 
if desired, only after all calculations are done. 

Since we are going to find the densities of measure of many different geo¬ 
metrical figures or bodies, we shall always use the symbol p for it, which will 
refer to different figures or bodies in different sections. Sometimes different 
densities are used in the same section. In that case either subscripts will be 
used or else they may be distinguish by their arguments. For example in sec¬ 
tion 4 p(A, A', B, B'\ f) is the value of the measure density of the length of the 
chords that join the segments AA! and BB' at £, whereas p(OA',OB', 7 ; f) is 
the value of the measure density of the length of the chords that join the seg¬ 
ments OA' and OB' (where ZA'OB 1 = 7 ) at i. Still to avoid heavy notation, 
given two points A and B, by AB we may mean either a segment of extremes 
A and B or the length of the said segment. 

Lemma 2.1. The longest chord parallel to any given direction hits one of 
the vertices of the triangle. 

Proof. Since the three vertices of a triangle are not aligned, as we move 
along a direction perpendicular to the set of parallel straight lines, we meet 
the three vertices at different times. When the first vertex is met the length of 
the chord is zero. Then it increases until the second vertex is met, after which 
it decreases again to be of length zero at the third vertex. □ 

Lemma 2.2. The length of the chords of a triangle parallel to any given 
direction is uniformly distributed between zero and its maximum length. 

Proof. Consider the set of chords of a triangle parallel to any given direc¬ 
tion. Their length changes linearly with the displacement as we move perpen¬ 
dicularly to that direction. □ 

Definition 2.1 We say that a chord of a given direction is “dominated” by 
the vertex A when chords having that direction reach their maximum length 
when they meet the vertex A. We denote this set by D(A). 

Lemma 2.3. The chords of a triangle belong to D(A)U D(B)U D(C). For 
any two vertices A and B, p(D(A) D D(B )) = 0. 


4 


Proof. The first statement is obvious. As for the second one, there is only 
one direction in D(A) D D(B). □ 

We now proceed to find the distribution of measure for the length of the 
chords dominated by a vertex, and later we may add the contributions of each 
vertex. 



Figure 1: Vertices A and B are both acute 


When computing the distribution of measure of D(C), vertices A and B 
may or may not be both acute. We consider these two cases separately. 

In the first case, the set D{C ) is the set of chords whose direction 9 e 
[a, 7 r — /3], as seen in the figure. We denote by D(C,9) the set of chords of 
direction 6 which are in D(C), and by pc(9) its measure (restricted to the 
direction 9). This measure is simply the cross section: 


Hc(8) — csind. 


We note that 


= 




d9 c sin 9 = c(cos a + cos j3), 


which, by Cauchy’s theorem [ 8 j [HJ [30] [38], yields the identity 

c(cos a + cos j3) + 6 (cos 7 + cos a ) + a(cos j3 + cos 7 ) = a + b + c, 


( 1 ) 

( 2 ) 

(3) 


which can be obtained by direct and simple means. 

We denote by he (9) the length of the longest chord in D(Cj9). 


h c {0) 


he 

sin 9 


(4) 


The measure density of the length of the chords in D(C, 9) is constant when 
the length i is smaller than hc{9 ), and zero otherwise. Since the measure of 
D(C,9 ) is ^c{9), the measure density (which we denote by pc) is 
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(5) 


p C (6,£) = ^H(h c (6)-e), 

where H is the Heaviside function, which is 1 when its argument is positive and 
0 otherwise. The measure density pc(9,Z) contributes only to those Z which 
are smaller than hc(9). This is seen in Figure El which depicts the measure 
densities for two different angles in [a,n — (3\ which have different he (9). Only 
the measure density pc(9 2 , Z) contributes to the measure density at the Z shown 
in Figure El which is between he(9 1 ) and hc(0 2 ). 

Formally, the distribution of measure for the length of the chords in D(C ) 
is 


hi ell hi 

Figure 2: Density of measure for the length of the chords in D(C,9 1 ) and 

D(c,e 2 ) 


pc(() = jT Oe - 1). (6) 

We are going to substitute the Heaviside function in the integrand by appropri¬ 
ate limits for the integral for each l. When l < he, all the pc{0, Z) contribute 
to pc(Z). When he < Z < min(a, b ), consider the two chords of length Z which 
meet at the vertex C. They form angles arcsin and n — arcsin with the 
side c, where by arcsin we mean the value between 0 and n/2. Then it is 
clear that the chords of length in [he, min(a, b )] belong only to the sets D(C, 9) 
s. t. 9 e [a, arcsin U [n — arcsin — /3\. Likewise, if, say, a < b, and 
a < Z < b, then the chords are in the sets D(C, 9) s. t. 9 e [a, arcsin ^f-\. Then 

(7) 
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To do the integral we use 


, 2A 
he — —, 

c 

where A is the area of the triangle, formulae ([Tj) and (HJ) which yield 

^ = A <*n*e 


h c (0) 2 A 


and the def ini te integral 

pbb 


sin" 9 = - ( bb — aa + 


sin 2 aa — sin 2 bb 


( 8 ) 


( 9 ) 


( 10 ) 


The result is 


Pc(i) = 


4 A 


(2A 




sin 2a sin 2/3 


+ 


+ 


TT/n 2A. tt . 1 n . ( .2 A sin 2a 2A I [2A\ . 

H(t - —)H(b - l) ( aresm - - a + — --J1 - l-\ I + 


9 A 

HU - )H(a 

c 


2A n sin 2/3 2 A I (2A\ 

aresm —-— pH---r a /1 — I —- I 


d 


d 


\d J 


• ( 11 ) 


Using the definition 


PaAA-, Oi, (3, 7 , a, b, c) = p c (Ph ( 12 ) 

we may write the measure density of length of chords of an acute triangle as 


Pac.tr(<) = Pzc((, P, 7.«. c, o) + p K (<, 7, a, ft c. a,b,) + p^((, a, ft, 7, o, 6, c). 

(13) 

We could, of course, have omitted the arguments a, (3 and 7 , because they are 
are determined by a, b and c, but we believe things are clearer this way. 

We now consider the case in which one of the vertices opposite to the vertex 
considered is obtuse. In Figure [3] the same obtuse triangle as in Figure |T] is 
considered, but now it is lying on side a. All reasonings leading to formula ([ 6 ]) 
still hold. This formula for vertex A is 
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Figure 3: Vertices B and C are not both acute 


Pa(() = f 1 it> ^j»(M») - 4 (14) 

The reasonings after formula (J 6 J) yield now only two terms, because the entire 
triangle lays now to only one side of the height h A . 


PA (e) = H(b-e ) 



Va(6) 

h A (0) 


+ H(i-b)H(c 





Va{0) 

h A (0) 


a 


2 


4 A 


H (b — £) 



+ 


sin 2 7 \ 
~ 2 ~ ) 


+ 


H(JL — b)H(c — £) 


( 2 A 

arcsm —- — 

l ai 


P + 


sin 2/3 
2 



(15) 


Using the definition 


PM(l,a,P,1,a,b,c) = PaW, 


(16) 


we may write the measure density of length of chords of an obtuse triangle as 


Poht.tM = Pm((,o, +Pobi((,P, a, 7, b, a, c) + p^(e,a,/3,j,a,b,c), 

(17) 


where c is the longest side. 














We have proven 

Proposition 2.1. The density of measure of the length of chords of a 
triangle, when the set chords is endowed with the measure p, is given by ex¬ 
pression when the triangle has an obtuse angle, and by expression (C3|) 
otherwise. 

Recently Gasparyan and Ohanyan |20] have found the density of measure 
of the length of chords of a triangle. They have not, however, organized their 
results in a form as compact as here, where all that is needed to find the 
density of measure of the length of chords of a triangle is formulae (fTTlh(fT3l) 
and (PSD-CETD. 


3 Distribution of the length of chords of two 
concurrent segments 

Let us denote by Cr(a ) the set of chords which cross the side a of a triangle. 
From this definition these three statements follow: 1) D(A) C Cr(a), D(B) C 
Cr(b ) and D(C ) C Cr(c); 2) up to sets of measure zero, each chord is contained 
in exactly two of the sets Cr(a ), Cr(b ) and C'r(c); 3) D(C) = (D(C)r\Cr(a)) U 
(D(C)nCr(b)),D(B) = (D{B)nCr(a))U(D(B)nCr(c)) and D{A) = (D{A )n 
Cr(b)) U ( D(A ) flCV(c)), where U denotes disjoint (up to sets of measure zero) 
union. The relations among the sets H(vertex) and CV(side), which were 
summarized at the beginning of these section, are depicted in Figure [H 


Cr(a) 



Figure 4: Relations among the sets H(vertex) and CV(side) 


We are going to find the distribution of the length of the chords of the sets 
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of the form D(C ) D Cr(a). For each direction 9 of the chords, the density 
of measure of D(C,9 ) fl Cr(a,9) and D(C,9 ) fl Cr(b,9) is a step function of 
the same width as the density of measure of D(C, 9 ), but with shorter heights 
which add up to the height of the density of measure of D(C,9). In Figure Q] 
we saw that this height was csin 9 (formula ([I])). In the same figure we can 
see that, in an obvious notation, 


(ji{D{C,6)nCr{b,6)) = bsm(9-a) (18) 

and 

n(D(C, 9) fl Cr(a, 9)) = asin( 7 r — /3 — 9). (19) 

We now follow the steps of section 1 changing pc{9) and Pa(9) by the 
above expressions. 

For D(C,9 ) fl Cr(b,9 ) in formula (5) we substitute pc{9) by //( D(C,9 ) fl 
Cr(b,9)), and in formulae (9) and (10) we substitute their right hand sides 

by || sin 9 sin(# — a) and | ^( bb — aa ) cos a + sin(aa — bb ) cos(aa + bb — a) j, 
respectively. We prefer to call the two segments which meet a and b, rather 
than b and c. Therefore we do the substitution (a, a) —> (c, 7 ), (6, /3) (a, a) 

and (c, 7 ) —» (b, (3) to find that the density of measure in Cr(a) fl D(B ) is 


= Ta 


H1 T 


£ ) [p cos 7 + cos a sin /3) + 


H{c-£)H\t- 


2 A 


2 A 2 A I (2Ay \ (2A\ . 

arcsm —-— yl— ( — ) —a ) cos 7 +( — ] sm 7 —sin a cos p 


be) ~) ' 1 \be 


H(a — e)H\ e 


2 A 


2 A 2 A , 

^Te-Te ' 11 


2 A 

~bl 


7 cos 7 +I 


2/1Y 


- 


be) 


( 20 ) 

When the height tig does not hit the side b, as shown in Figure 0 there 
are only two terms: 

ab 

PaB (e) = — [H (c - £) ({3 cos 7 + cos a sin /?) + 


sm 


HV-oma-ti 11 — ^-^V 1 -(^) 2 - T ) C03T+ ( 1 -(^) 2|3m1 ' 


When 7 > f (Figure EJ), a similar calculation yields 
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B 



Figure 5: The height Jib does not hit the side b 


B 



Figure 6 : 7 > | 


PaB{£) — 


ab 

2A 


H(a -£) 


/■*-“ ah 

I dd sindsin(d— 7 )+—.£/"(£—a)id(c—£) 


' 7T—arcsm 


^A 


dd sindsin(d— 7 ) = 


ab 

4A 


H (a — £) (/3 cos 7 + cos a sin /?) + 


H(l—a)H(c—£) 


2 A 


arcsm — — 
bi 


2 A 

n 




+sina! cos a —a 


cos 7 + 



2 


(sin a ) 


sm 7 


Since 


Cr(a ) n Cr{b) = (Cr(a) D D{B)) U (Cr{b) D D(A)), (21) 
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to find the density of measure of the length of the chords which cross the sides 
a and b we have to add the corresponding densities of the chords in the sets 
Cr(a) n D(B) and Cr{b) D D(A). 

We have proven 

Proposition 3.1. The density of measure of the length of chords which 
cut two concurring segments of lengths a and b and which form an angle 7 is 
given by 


p(a, b, 7 ; £) = paBif, a, P, 7 , a , b, c) + p bA {£, a, P, 7 , a, b , c ). ( 22 ) 

4 Distribution of the length of chords between 
two segments 

Proposition 4.1. The density of measure of the length of chords which cut 
two non parallel segments of extremes A, A' and B, B' and which form an angle 
7 7^ 0 is given at £ by 


p{A, A 1 , B, B'] £) — 

p{OA', OB' , 7; £)- plOA, OB' , 7; £)-p{OA!, OB , 7; £) +p(OA, OB, 7; £). ( 23 ) 

where the p’s in the rhs denote the density of measure of Proposition 3.1 and 
OA, OA’, OB, OB’ denote lengths. 



gamma a__A’ 

Figure 7 : Any two non parallel segments define an angle 


Proof. It is clear from Figure |T| that 

plOAf, OB', 7; t) = p(OA, OB', 7; £)+p(OA', OB, 7; £)-p(OA, OB, 7; £)+p(A, A', B, B'- £), 

from which the result follows immediately. □ 

The parallel segments case can be found as a limit of the previous case, but 
it is faster to obtain it directly, because when the segments are parallel there 
are, for a given length, at most two directions for which the chords have that 
given length. 
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B’ 


v theta 


s pi - theta 


A’ 


Figure 8: Two parallel segments. 


In Figure [8] we see two segments whose supporting lines are a distance 
d apart. The density of measure of the chords of length £ = d is obviously 
proportional to the length of the intersection between AA' and the projection 
of BB' onto the x axis, that is, 


where 


p(d) oc h(y mines', xa') — max^,!^)), 
h(x) = xH(x). 


(24) 

(25) 


When £ > d there are two directions for which the chords have length £. These 
directions are arcsin and 7r — arcsin The density of measure of the chords 
of length £ and direction arcsin j (denoted 6 in Figure |H|) is proportional to the 
length of the intersection between AA' and the projection of BB' along the 
lines of direction arcsin j times j. A factor j is needed to convert the length 
of the intersection into the cross section of the set of chords. Finally, an angle 
dependent factor is needed to compute the density of measure. Since 


6 = arcsin 


d 


dd = 


d 


zd£, 


(26) 


this factor is 


d 


£^/£ 2 - d 2 

The measure of the chords of length l and direction 


7r — arcsin | can be computed similarly. 


pW 


d 2 

£ 2 V £ 2 - d 2 


h(mm(xB' 


V £ 2 — d 2 , xa ') — max(i B — \J £ 2 — d 2 , x^)) + 


h( min (xb> + V£ 2 — d 2 , xa i ) — ma x(xb + \/£ 2 — d 2 , xa )) 


> d. (27) 
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If the previous definition is taken to hold also when £ = d, only a measure zero 
error is made. 

We have proven 

Proposition 4.2. The density of measure of the length of chords which 
cut two parallel segments whose supporting lines are at a distance d is given 
by the previous expression when £ > d and is zero otherwise. 


5 Distribution of chords in a convex polygon 

In order to compute the first formula of the preceding section when A, A', B 
and B' are given in Cartesian coordinates we need to compute, first of all, the 
coordinates of O. It is tedious but elementary to find that the straight lines 
supporting the segments AA 1 and BB' meet at a point O of coordinates 


(x 0 ,lJo) 


XA 

VA 

XA' 

VA' 

X B 

Vb 

X B ' 

Vb' 

XA' 

- X A 

X B ' 

- x B 


xa> ~ X A 


xb> ~ x B 


VA’ ~ IJA 
UB> ~ IJB 


XA 

VA 

XA' 

VA' 

x B 

Vb 

x B > 

IJB' 


XA> - X A 
Xb< - X B 


VA ' - VA 


Ub’ ~ Vb 


VA’ ~ VA 
Ub> ~ Vb 


(28) 


V 


/ 


The lengths which appear in formulae (22) (to which formula (23) refers) 
and (23) are OA, OA', OB, OB' and A'B', all of which are straightforward to 
compute because the coordinates of the five points involved are now known. 
As for the angles of the triangle OA'B' , they are 


AA' ■ BB' 

7 = arccos ) -p, 

\AA!\\BB'\ 


and, by the cosine theorem, 


a = arccos 


\BB'\ - \AA!\ cosy 


|kL4'| 2 + \BB'\ 2 - 2\AA'\\BB'\ cosy 


(29) 


and 
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( 32 ) 


_ \AA'\ — \BB'\cosp 

p = arccos — , = 

^\Al'\ 2 + \BB'\ 2 - 2\Al'\\BB , \cosp 

The “missing” side of the triangle can also be found using the cosine theorem. 

It is now possible to write p(A,A',B,B';£) in Cartesian coordinates. 

Proposition 5.1. The density of measure of the length of chords of a 
convex polygon given by the Cartesian coordinates of its vertices can be found 
using the formulae of sections 4 and 5. 

Proof. In a convex polygon, up to measure zero, all chords go trough 
exactly two sides. If a convex polygon is given by the Cartesian coordinates 
of its vertices, then its total measure density is the sum of the densities of all 
pairs of sides, which the formulae of this and the preceding section show how 
to find. □ 

That the set of chords cutting a convex polygon can be decomposed in the 
(up to measure zero) disjoint union of the sets of chords cutting the pairs of 
sides is a straightforward proposition which has been known and used by other 
authors (e. g. [33, [29]). The advantage of the formulae presented up to here 
over the formulae presented by H. S. Harutyunyan and V. K. Ohanyanian [29] 
is that the formulae of sections 4 and 5 allow for a straightforward substitution 
of the Cartesian coordinates (cf. [29], section 5). The author has written a 
file of Mathematica code that does just that and which can be downloaded at 

H2I- 

6 Convex examples 

We have shown how to compute the distribution of chords of a (convex) poly¬ 
gon by computing it for each pair and then adding all the distributions found 
in this manner. This leads naturally to the question of whether a couple of 
non-congruent n-gons exists such that the (") pairs of sides are, setwise, con¬ 
gruent to each other. Such a pair (actually, a whole family of pairs) was found 
by Mallows and Clark [32] in 1970, and we show in Figure [9] their measure 
density. 

The particular pair taken as example is based on the octagon of coordinates 

{A,B,C,D, 

E,F,G,H} = {( 1,0), (1/V2,1/V2), (0,1), (-1/^2,1/^2), (-1,0), (-1/^2,-1/^2), 
(0, —1), (l/\/2, — 1 /\/2)}. To obtain one of the dodecagons we add the ver¬ 
tices {AE>, BC, 

DE,EF} = {l.l^y^, 1.1^4^Z, 1.1^242^ l.l^yE}, and to obtain the other the 
vertices {AB, DE,EF, FG} = {1.1^, 1.1^±®, l.l!±£, l.iM}. Although 
not clear in Figure [91 the graph is continuous, since the formulae obtained 
in the preceding section are continuous. More generally, this follows from a 
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Figure 9: Measure density of a pair of Mallows and Clark dodecagons 

theorem of Sulanke for bounded, non-degenerate, convex figures ( 05 . page 55, 
Proposition 2). This particular example is worked out in the hie Convex.nb 
in [19} . With that code the above graph is produced in about 30 seconds on a 
personal computer. 

7 The concave quadrilateral 

The measure density for the depicted quadrilateral can be decomposed in mea¬ 
sure densities of pairs of segments as follows: 

P(V A,AO) T P(V D,DO) T P(VA,OC) T P(VD,OB) T P(AO,OC) T P(BO,OD) ~P(OB,OC) +Pconc(AB,CD) ■ 


(33) 


Of these eight terms the first seven are of the kind already investigated, while 
the last one is new, because not every segment joining AB and CD is contained 
in the quadrilateral. 

The new term can be approximated (and bounded by above) by 



p(M 0 , M 0 +rriiU, N 0 , op (M 0 + mjw); £)+p(M 0 +m 1 u, M 0 +m 2 u, N 0 , op (M 0 + m 2 u ); £) 
+p(M 0 +m 2 u, M 0 +m 3 u , N 0 , op (M 0 + m 3 u ); £)+p(M 0 +m 3 u , M 0 +m 4 u, N 0 , op (M 0 + m 4 u ); £), 


(34) 
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Figure 10: The concave quadrilateral. 


where u = , rn l = \M t — M 0 |, op is a geometrical function which assigns 

to a point in the left side the point in the right side such that the straight 
line which they determine contains the vertex V and p is the measure density 
function for two segments which has already been studied. The function op is 
derived in the Appendix. 

For any n e N there is a set of equally spaced points {M 0 , Ml, • • • , M n } in 
the right side and an approximation as above can be written: 


Y p(M 0 +- - -\M 1 -M 0 \u,M 0 +-\M 1 -M 0 \u,N 0 ,op (M 0 + -| M 1 - M 0 \u ) ;f). 

' n n \ n J 

i=i x x 

(35) 

This expression is the starting point for two different ways of computing 
the measure density for a concave quadrilateral. The first way to be presented 
is numerically more efficient, while the second one is mathematically more 
sophisticated and allows, sometimes, an analytical computation of the measure 
density. 

In section 10 we shall refer repeatedly to a pentagon of the kind ABCDV. 
For this reason we give the following definition. 

Definition 7.1 A concave pentagon whose concave vertex is the intersec¬ 
tion of the diagonals of the convex quadrilateral whose vertices are the other 
four points will be called a diagonal pentagon. 
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Note that there are four diagonal pentagons for a given set of five vertices 
satisfying the above definition. Note also that in a diagonal pentagon the 
distribution of the length of chords between any pair of sides can be computed 
using Propositions 4.1 and 4.2, except for one pair, which is the pair AB, CD 
in the case of Figure fTTTl 

8 Converging bounds for the measure density 
of the concave quadrilateral 

As remarked in the last section, the last expression is an upper bound for the 
measure density p CO nc(AB,CD)- This is better seen in the case n = 4, depicted 
in Figure (TUJ In the expression (153]) all chords joining the segments AB and 
CD and passing above the vertex are counted in, as well as some chords which 
pass below the vertex V. Similarly, the expression 

^ " 1 " f " 1 \ 

VVm 0 +-- \Mi—M 0 \u, Mq+—\Mi~Mo\u, No, op ( M 0 + -- \M X - M 0 \u ) -,£). 

n n \ n J 

2—1 x 7 

(36) 

is, for any n G N, a lower bound for the measure density p CO nc{AB,CD)- 

Using the notation app_(n)(£) and app + (n)(£) for the last two sums, we 
have remarked that 

app_(n)(£) < p conc (AB,CD )(£) < app + (n)(£), Vn 6 N. (37) 

The difference of the sums is 

app + (n)(£) — app_{n)(t) = 

n . _ 

y p(Mo H-1Mi — Mq\u, Ma H—|Mi — Mq\u, 

n n 

2=1 

op[M 0 + - — -\M X — M 0 \u) , op (M 0 +—\M X — M 0 \u) ; £) —» 0. (38) 

\ n J \ n J n—>oo 

This is intuitive because the difference app + (n) (£) — app_ (n) {£') involves, as 
n —> oo, only the chords that go through the vertex V, which is a set of 
measure zero. It can be checked rigorously using the analytic form of the 
density of measure for the chords between two segments. 

Proposition 8.1. The last two equations prove that 

lim app- (n) (£) = lim app + (n) (£) = p conc (AB,CD)(£). (39) 


18 






9 The measure density of the concave quadri¬ 
lateral as a Riemann integral 

We momentarily omit the arguments No and l of p in the sum (35) and use 
the following short hand notation for the other three: 

n 

^2p(i-l,i,op(i)). (40) 

%— 1 

p is C 2 except at a finite number of points (see formulae in the sections on 
convex polygons), therefore: 


n 

Y p ^ - 1 ,i,op(i)) = 

2—1 


71 \ B _ M / 1 \ 

Y p(i - 1, * - 1, op (*)) + p' 2 (i - 1, i - 1, op (i)) - -- + o ( - ) (41) 

n \n1 

i= x v 7 

where p 2 is the partial derivative with respect to the second argument. The 
term p(i — 1, i — 1, op (*)) is zero because the segment of extremes Mq + — \B — 
A\u and Mq + — | B — A\u is of zero length. The expression can be further 
transformed as follows: 


Y P'iii ~ 1, * - 1, op (i)) ——— + o 

Z — J n \n) 

2—1 X / 

Y P'iii ~l,i~l,op(i- 1)) ——— + o (- 

n \n 

i =1 v 

Pxsd -l,i-l,op(i- 1 ))op\i - 1) ^ + o 

n- 



YP'iii ^l,i-l,op(i- 1)) ——— + o . (42) 

Z ' n \ n J 

2— 1 X ' 

Had we started from expression ()36l) we would have also reached expression 
(42). Since expression (j36|) is, for any n, always smaller than or equal to 
expression (1361) and the differ by at most Y^i=i 0 (y), which satishes 



lim n o 

n—>• oo 



lo(i) 
lim n - 

n—y oo Tl — 


n 


lim 

n—>-oo 



n 


0, (43) 
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then the limit of expression (42) as n —> oo is the following Riemann integral: 


\B-A\ 


da p' 2 (a, a, op(a )) = 


\B-A\ 


da p' 2 {A + au, A + au , C, op (A + au ); £), 


(44) 

where in the rhs we drop the short hand notation and A, B and C are the 
extremes of the thick intervals depicted in Figure [TO] 

We have proven 

Proposition 9.1. 


\B-A\ 

da p' 2 (A + au, A + au, C, op (A + au ); £) = p conc (AB, CD){£). (45) 

In the remaining part of the section a formula for the integrand of the 
above expression will be derived. Note that the approximations app_[n)(£) 
and app + (n)(£) discussed in the previous section are lower and upper bounds 
of p 2 (£), but not lower and upper Riemann integrals of the preceding Riemann 
integral. 



9.1 The derivative p' 2 

We need to find the derivative p' 2 . In this subsection we use the notation of 
section 3. First we address the general case in which the two segments are not 
parallel. From (23) it is clear that 

P 2 (A, A', B, B 1 ; l) = p\(OA\ OB', 7 ; £) - p'^OA', OB, 7 ; £). (46) 

Since 

p(a, b, 7 ; £) = p aB (£, a, f3, 7 , a, b, c) + p bA (£, a, ft, 7 , a, b, c ), (22) 

we need to find the total derivatives Ap aB {£ , a, j3, 7 , a, b, c ) and Ap bA {£, a, (3 , 7 , a, b, c ). 
The first of these total derivatives is 

d_ 
da 

PaBbi^ A 7, a, b, c)da + p' aB7 {£, a, /3 , 7 , a, b, c)A-da 

From the cosine theorem, c = \Ja 2 + b 2 — 2 ab cos 7 , 

dc a — b cos 7 
da c 


(47) 

(48) 


dot d/3 

p aB (£, a, ft, 7 , a, b, c) = f> aB2 {£, a, ft, 7 , a, b, c)—da+p' aB3 (i, a, j3, 7 , a, b, c)— da+ 
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From the sine theorem, 

’ a 


Likewise, 


sin 7 


yj a?+b 2 — 2 ab cos 7 

da b sin 7 
da c 2 


a = arcsm a 


sin 7 


yj a?-\-b 2 — 2 ab cos 7 5 


(49) 


d/3 _ b sin 7 

da c 2 ' ^ 

A similar calculation is required for a, /3, 7 , a, b, c). This approach 

leads to an expression for p ' 2 which is a few pages long and, therefore, of 
questionable advantage with respect to simulations. 


b 



Figure 11: In the example above (3 ~ 7 t/2 . x rnax {f3 r ) is, for a given /3', the 
maximum x for which the chord still cuts the upper side. 


9.2 A more direct calculation of p' 2 . Preliminaries 

p' 2 can be found by a more direct method, which is to consider the distribution 
of the length of chords between two segments, one of which is of infinitesimal 
length. 

As a first step towards this goal we are going to find the measure of chords 
that go through two concurring segments as a double integral over the seg¬ 
ments. Let x and y be coordinates along the said segments, as shown in 
Figure [TTJ To find the said measure we find the section of the chords for a 
given direction and then integrate over the directions: 
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h = 



d/3' s((3') 




ds((3', x ) 
dx 



ds((3', x) 
dx 


d/3 1 / dx sin f3' + 


r» 7T —7 


d/3' 


. 09 ') 


dx sin/S', 


( 51 ) 


Jo Jo J p Jo 

where (3 is the angle between the horizontal side and the missing side of the 
triangle and 


ds = sin /3'dx, (52) 

which is obvious from Figure [TU has been used. 



Figure 12: For a given x, there is a maximum /3' max {x) for which the chord still 
cuts the upper side. 


To change the integration variable in the last integral from /3' to y we use 
the sine theorem in the form 


sin j3' 
V 


sin 7 

\Jx 1 + y 2 — 2 xy cos 7 


/3' = arcsin 


y sm 7 


\Jx 1 + y 2 — 2 xy cos 7 


d/3' x sin 7 

dy x 2 + y 2 — 2 xy cos 7 


(54) 
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Substitution in the last integral yields 


P = 



dx dy 


xy( sin 7 )" 


0 Jo 


(. x 2 + y 2 — 2 xy cos y ) 3 / 2 


(55) 


The above formula holds for two concurring segments. 

To take the 7 —> 0 limit in the above integrand we note that the denom¬ 
inator is the cube of the distance between the points of coordinates x and y 
and that the numerator is the product of the heights of the triangle which go 
through the points of coordinates x and y. Then, when the lines are parallel, 
the above integrand becomes 


(d 2 + (y — x ) 2 ) 3 / 2 ’ v ; 

where d is the distance between the parallel lines. 

9.3 A more direct calculation of p' 2 . Conclusion 

It follows from the equality (1451) that p' 2 (A + aft, A + aft, C, op {A + aft); £) da 
is the measure density of l for the chords which join the segments A + aft, A + 
(a + da)u and C, op ( A + au ). Therefore, using a notation that keeps only the 
last argument of p' 2 , we may compute it as follows: 



Figure 13: 


m 


lim — 

dx-y 0 dx 




xy( siny ) 2 

(. x 2 + y 2 — 2 xy cos y) 3 / 2 


d(\/x 2 + y 2 


2 xy cosy— i), 


. (57) 

where the coordinates x and y have their origin in the vertex shown in Figure 
IT3l The comparison of Figures 10 and 13 yields the relations between the 
arguments of p' 2 in (1451) and the arguments of p 2 used above: 
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\0C\ —> yi, |Oop (A + au) | —» ?/ 2 , |0A — em| —> x. (58) 

To take advantage of the Dirac delta function we have to change the vari¬ 
ables of integration from (x, y) to <p(x, y) = y/x 2 + y 2 — 2 xy cos 7 and some 
other variable. We choose the other one to be the angle /3'. By the sine 
theorem, the relation between the variables is: 


x = sin (7 + j3') 

sin 7 \ 1 ~ / 


y 

The integrand becomes 


-4— sin fj' 

sin 7 


(59) 


d(x,y) 

d (4,4') 


xy( siny ) 2 

(. x 2 + y 2 — 2 x|/cos 7) 3 / 2 


5( a/x 2 + y 2 


2 xy cos 7 — £) = 


cf) sin (3 r sin( / 9 / + 7) 
sin 7 ( j ) 


8(4-1) 


sin f3' sin( / 5 / + 7) 
sin 7 


6(4-e). 


(60) 


We are going to integrate first over <fi. For hxed /3' the limits of integration of 
(j) are, by the sine theorem, sin s ^,^ 7 ^ x and sin s ^,^ 7 ^ (x + dx). Therefore, in the 
expression (where the integration limits have been omitted) 


m 


lim — 

dx^f 0 dx 


0 sin/3 / sin(/3 / + 7) 
sin 7 


d(j) 8(4 — £), 


(61) 


the integral over 4 will be 1 when £ G 


sin 7 ^ 

sin(/?'+7) 


sin 7 

sin(/3'+7) 


(x + dx) 


and 


0 otherwise. The way to capture this into the above expression is to give 
appropriate limits of integration in the integral over fd'. To do this we note 


that 


sm 7 sm 7 

x, ^ (x + dx) 


Lsin( / 9 / + 7) ’ sin( / 9 / + 7) 


(62) 


is equivalent to 


/3'e 


X 

arcsm I — sin 7 I — 7, arcsm 


x + dx 

—-—S1117 1 -7 


U 


7 r — arcsin 


x + dx 


sm 7 


— 7, 7 r — arcsin 


x 


sm 7 



h U I 2 , (63) 


where I\ and I 2 are notations for the intervals above and by arcsin the branch 
that takes values between —7 t/2 and 7 t/2 is meant. The two intervals corre¬ 
spond to the two chords that have a given length £, one to each side of the 
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shortest segment that goes from the point of coordinate x to the straight line 
where the point of coordinate y lies. Let the angles (3' which correspond to 
the coordinates y\ and ?/ 2 be and /3 2 (which satisfy sin /hi = y sin 7 and 
sin fj 2 = y sin 7 , respectively). Therefore, 


m 


lim -L f d 0 + + llm _L f ¥ sjnlMl+j} 

dx^o dx J hn[ p lM sm 7 dx-+o dx J hn[l3lM sm 7 


(64) 

Since the intervals I\ and J 2 are infinitesimal, they can be made so small 
that they either are fully contained in [/hi, /h 2 ] or are disjoint with it. Sup¬ 
pose that 1 1 C [/3i, /3 2 ]. Then the hrst integral is, to hrst order in dx , the 
value of its integrand when /3' = arcsin (| sin 7 ) — 7 times the width of 
the interval I \. When f3' = arcsin sin 7 ) — 7 , then sin (3' sin( / S / + 7 ) = 
sin (arcsin (|sin 7 ) — 7 ) |sin 7 . To find the width of the interval I\ we use 
the expansion 


x dx \ / x 

arcsm | --— sin 71—7 = arcsin ( — sin 7 ) — 7 + 


S1117 


\Jt 2 — x 2 sin 2 7 


Therefore, if I\ C [/3i,/ 3 2 ], the first term of expression (64) is 


x sm 7 


X 

sm I arcsm I — sin 7 1 — 7 


i\J(? — x 2 sin 2 7 
Likewise, if J 2 C [/3i,/ 3 2 ], the second term of expression (64) is 


x sm 7 


t\/1 2 — x 2 sin 2 7 


x 


sm I arcsm — sm 7 1 + 7 


dx + o(dx). 

(65) 

( 66 ) 

(67) 


The sine theorem yields y = . The y coordinates which correspond 

to /T = arcsin (| sin 7 ) — 7 and to (3' = n — arcsin (| sin 7 ) — 7 are 


and 


y- = 


y+ = 


sin ( arcsin I f sm 71—7) 


sm7 


in arcsin (j sin 
sin 7 


( 68 ) 


respectively (note that the functions y± take a complex value when x sin 7 > £, 
i. e., when t is too short to reach the y side). Therefore, 


P 2 W — 1 [l/l , 2 / 2 ] ( 2 /— y 


X sin 7 


£\/£ 2 — x 2 sin 2 7 


x 


sm I arcsm I — sm 7 ) — 7 ) + 
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^[yi,V2](y+) 


x sm 7 


l\Jl 2 — x 2 sin 2 7 


x 


sm I arcsm I — sm 7 I + 7 


where is the indicator function of the set [2/1,2/2], that is, 1 [yi,y 2 ]{y) 

if y G [ 7 / 1 , 7 / 2 ] and 0 otherwise. 


(69) 
= 1 


y2.-'' 



Figure 14: In this case 1 [ yi , y2 \(y-) = 0 and l[ yi , y2 ](y+) = 1 


When 7 < 7t/2 there are two segments of length i that go from the point of 
coordinate x to the straight line supporting the coordinate y. When 7 > 7t/2 
there is only one such segment. Consequently, when 7 < 7t/2, each of the 
prefactors may be 0 or 1, while if 7 > 7t/2, the first prefactor is always 0. 

When 7 = 0 and the distance between the parallel lines is d, the remarks 
previous to formula (156|) lead to 


and 


t/_ = x — \/i 2 — d 2 , 
y + = x + \Jt 2 — d 2 , 


(70) 


f ( \ cP 

P2W — -h [7/1,2/2] (7/— ) + l[j/i,y 2 ] (y+) J 


d- 2 


(71) 


To use the preceding four formulae to compute the density of measure as 
prescribed in formula fl45l) . we note that the angle 7 is the angle AOD of Figure 
HOl and we reverse the change in notation (l58|l : 


7/1 —> \OC\, 7/2 —> \Oop {A + au) |, x—*\OA — au\. (72) 

We have proven 

Proposition 9.2. The integrand p' 2 of Proposition 9.1 is given by formulae 
(ESP and If69 ]). The angle 7 is the angle AOD of Figure ITdl and the change of 
notation [58 1) is reversed: 


7/1 —> |OC|, 7/2 —> |Oop(A + au) I, x—}\OA — au\. (73) 
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When 7 = 0 and the distance between the parallel lines is d, then the integrand 
p' 2 of Proposition 9.1 is given by formulae HM) and 


10 The general case 

The quadrilateral in section 7 had four sides to start with. But in order to 
compute its cld, its sides OA and OD were conveniently divided into two seg¬ 
ments. Together with the sides VA and VD this gave a total of six segments, 
which gave 15 pairs. Five of the sets of chords joining the segments of each 
pair were empty, nine were of the kind studied in sections 2-6, and the last set 
of chords was of a new kind whose cld was derived in sections 7-9. Thus the 
non empty sets of chords between two of these segments are either the chords 
between two opposing sides of a convex quadrilateral or the chords between two 
particular sides of a diagonal pentagon. In this section we prove that such a 
procedure can be applied to any polygon, that is, the concavity can always be 
packed in diagonal pentagons, so to speak. 

Definitions. Let AB and CD be two segments such that their convex cover 
is a quadrilateral. Then by Qu(AB , CD) we denote the set of chords joining 
the sides AB and CD of the said convex quadrilateral. Let AB and CD be 
two segments and V be a point such that BAVDC is a diagonal pentagon 
whose concave vertex is V. Then by DPent(BA, AVD, DC) we denote the set 
of chords joining the sides AB and CD of the diagonal pentagon BAVDC. 
When AB and CD are sides of a polygon, the segments in Qu(AB , CD) may 
or may not be fully contained in the polygon. The latter segments constitute 
the set of chords Ch(AB,CD). By definition, Ch(AB,CD) Q Qu(AB,CD). 

When the polygon is concave there is at least one pair of sides AB and CD 
for which 0 C Ch(AB,CD ) C Qu(AB,CD) holds, where the inclusions are 
not trivial. Then there is at least one segment s e Qu(AB , CD) which is a 
chord, s divides the convex quadrilateral ABCD into two parts (we suppose 
here that ABCD is the order in which the vertices yield a convex quadrilateral). 
At least one of the parts must contain a point P which is not in the polygon. 
Any segment joining the chord s and the point P must cross the perimeter of 
the polygon at least once. Therefore, in the said half of the quadrilateral lays 
a portion of the perimeter of the polygon. This portion cannot be just a part 
of a side, for in that case it would cross at least one of the sides AB and CD, 
which would not be sides themselves. Thus it must contain at least one vertex. 

Definitions. Consider the set which is the union of the points {A, D} and 
the portion of the perimeter of the polygon which is contained in the quadri¬ 
lateral ABCD and crosses the segment AD. We denote its convex cover by 
Per(A, D). Per(B,C) is defined likewise. By AB V we denote the intersection 
of the side AB and the chords in Ch(AB,CD). Informally speaking, it is the 
part of AB from which some point of CD is visible. CD V is defined likewise. 
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Proposition 10.1 For any two sides AB and CD of a polygon, the side 
AB can be covered by a finite collection of non overlapping segments with 
the following property. Let PaPb be any of the said segments. Then the set of 
chords that join PaPb and CD is either empty or is the set union or difference 
of at most three sets of chords of the form Qu or DPent defined above. 

Proof. Take any point P in the interior of AB V . Consider the chords 
Si(P) and s 2 (P) which meet the point P and are tangent to Per(A,D) and 
Per(B,C), respectively. Almost surely Si(P) and s 2 (P) touch only one vertex 
V\ of Per (A, D) and one vertex V 2 of Per (B,C), respectively. These are the 
vertices which limit the vision of the side CD from the point P. As the point 
P moves to either side it reaches points Pa and Pb such that the set of vertices 
touched by s/P) and s 2 (P ) increases by one vertex, because s/P) or s 2 (P) 
include some side of Per [A, D) or Per (B,C). Outside the segment PaPb the 
vision of the side CD is limited by a different pair of vertices. One can proceed 
in this manner and cover AB V by a set of disjoint segments such that in each 
of them the vision of the side CD is limited by a different pair of vertices. 

For each such segment PaPb, there are eight possible orderings of the points 
Qa, Q'a, Qb, Q'b> and we shall show that for each of them Proposition 10.1 is 
true. 

Case 1. In the case shown in Figure \TR Ch(P A PB,QAQ' B ) — 

Ch(P A P B , QaQb) + Ch{P A P B , Q b Q' a ) + Ch[P A P B , Q'aQ'b) = 

DPent(P A P B , PbV 2 Qa, QaQb)+Qu(PaPb, QbQ^PDP ent(P B PA, PaViQ'bi Q'bQ'a)- 

When Qa = Qb and/or Q' A = Q’ B two/one of the three previous terms 
disappear and the Proposition 10.1 holds. 



Figure 15 : Case 1 of Proposition 10.1 


Case 2. In the case shown in Figure 1771 C^PaPbjQaQ'b) = 

Ch(P A P B , QaQb) + CI/PaPb , Q'aQ'b) - Ch(P A P B , Q'aQb) = 

{DPcu/PaPb, PbV 2 Qa,QaQb)~ segments that join PaPb and QaQb and 
go below the vertex Vi) + (DPent(PsPA, PaV\Q' b , Q'bQ'a)~ se g men ts that join 









Figure 16: Case 2 of Proposition 10.1 


PaPb and Q' A Q' B and go above the vertex V 2 ) — (Qu(PaPb, Q'aQb) — segments 
that join P a Pb and Q' a Qb and go below the vertex V\—segments that join 
PaPb and Q'aQb and go above the vertex V 2 ). 

But the segments that join PaPb and QaQb and go below the vertex Vj 
are the segments that join P a Pb and Q' a Qb and go below the vertex Vj, and 
the segments that join P a Pb and Q' a Q'b an> d 9° above the vertex V 2 are the 
segments that join P a Pb and Q' a Qb and go above the vertex V 2 . Therefore, 
in the preceding equality between sets, the sets that have been named using 
sentences in English cancel and we obtain 

Ch(P A P B ,QAQ ' B ) = 

Ch(P A P B , QaQb ) + Ch{P A P B , Q'aQ'b ) - Ch{P A P B , Q'aQb) = 

DPent(PAPB 5 PbV 2 Qa, QaQb)+DP ent(P B P A , P A ViQ ' b , Q' b Q' a )—Qu(PaPb, Q'aQb)- 

Case 3. Q' A = Qb- This is a limiting case of both case 1 and case 2. Either 
way we obtain: 

Ch(P A PB, QaQ' b ) = DPent(P A P B , PbV 2 Qa, QAQB)+DPent(P B P A , PaViQ'b , Q'bQ'a)- 

Case 4- Q'a = Qa- In this case P A ,Vi, V 2 and Q A = Q' A are aligned. 
Ch(P A P B ,Q A Q' B ) decomposes as in Case 2. 

Case 5. Q' B = Q B . In this case P B ,V 2 ,V\ and Q B = Q' B are aligned. 

Ch(P A P B , Q a Q' b ) decomposes as in Case 2. 

Case 6. Q A = Q B . In this case V 2 is on the side CD. Ch(P A P B , Q A Q' B ) 
decomposes as in Case 1 but has only two terms. 

Case 7. Q' A = Q' B . In this case Vj is on the side CD. Ch(P A P B , Q A Q' B ) 
decomposes as in Case 1 but has only two terms. 

Case 8. Q A = Qb and Q' A = Q' B . In this case Vj and V 2 are on the side 
CD. Ch(P A P B , QaQ's) decomposes as in Case 1 but has only the Qu(P A P B , QbQ' a ) 
term. 

When the polygon is not simply connected, there might be holes in the in¬ 
terior of ABCD. Each of this holes has a convex cover which we denote by 
Per(Hole). Then the vision of the side CD from a point P e AB may consist 
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of more than one segment (up to h + 1 segments if h holes are caught in the 
quadrilateral ABCD) and instead of two chords Si{P) and s 2 (-P) there might 
now be up to 2 (h + 1) chords which are tangent to two vertices in Per(A, D ) 
and Per(B,C) and to 2 h vertices in Per(Hole). But each of the segments 
QaQ'b corresponding to each pair of vertices can still only fall in any of the 
eight cases considered above. Thus the proposition also holds for not simply 
connected polygons. 

Since the number of vertices in Per(A, D ) or Per(B,C) or Per(Hole) is 
finite, the number of segments PaPb is finite. □ 

11 Concave examples 

11.1 Concave pentagons 



Figure 17: A concave symmetric pentagon 


The measure distribution of the shown symmetric diagonal pentagon (Fig¬ 


ure 


77\) is 


2 PAB,BC + 2 pAB,CD + 2 PBC,CD + pconc BC,ED ~ 


2p convex AB.BC A 2 p convex AB,CD A 2 pconvex BC,CD A ^Pconvex BC,EDi (74) 

because for each chord c joining the sides BC and ED contained in the pen¬ 
tagon there is a mirror image of it, c!, joining the sides BC and ED but not 
contained in the pentagon. 
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To compute p conc bc,ed we may also use (70) and (71). If we let B,C,D 
and E be the vertices of a square of side length 1, then d = 1, y i = 0 and 
y- 2 (x) = l—x. Then (70) and (71) become 


y- — x — yJT 1 — 1, 

y+ = x + \/l 2 - 1 , 

and 

f>2^) = (l[o,i-x](2/-) + l[o,i-x](z/+)) 

Thus 


(75) 


(76) 


Pconc BC,Ed(?) = j dx (l [0A _ x] (x-Vl 2 ~ l) + l[0,l- 8 ] (x+V £ 2 - 1 )) - 


i (i . . , ^ 1 1 - \fl 2 ~ 1 

dx ^l[VF3i :(1+ vF3T) /2] (a;) + 1 [0 ) (i-v7^i)/2 ] Wj ~ py/ p = f ' ’ 

(77) 

w/ien £ G [1, \/2] and 0 otherwise. 

For every chord c joining BC and ED and passing below the vertex A, 
there is a symmetrical chord c! of the same length passing above the vertex A. 
Therefore using p conc bc,ed = \pconvcx bc,ed one obtains the same result as 
above. 

Another example is the concave non symmetric, non diagonal pentagon 
depicted in Figure Its measure density has terms which are computed as 
explained in the part of the article devoted to convex polygons: 



Pab,bc+Pab,cb'+Pbc,cb'+Pd'c,b'd+Pe'c,de+Pe'c,ea+Pcd,de+Pcd,ea+Pde,eA: 
and the following two terms, which are computed using (68) and (69): 

Pconc BD',B'D + Pconc D'E',DE 

When computing p conc bd’.b’d, 7 = 7r/2 . This implies that y± = ±V ‘ I 2 — x 2 
and that the first term of formula (69) vanishes (as remarked after formula 
(69)), leaving 


X 

P2W = \y\,yi\~jjfi ■ 


(78) 


31 




















Figure 18: A concave non symmetric pentagon 


If we let the coordinate x be the one along the segment CD, with origin in C, 
then 


r\CD\ _ x 

Pconc BD',B'd(£) = dx 1 [CD',Cop(x)] (79) 

J\CB'\ I 

In this case the indicator function is 1 when \CD'\ < \JI 2 — x 2 < , where 

xa and ua are the Cartesian coordinates of A with respect to a system with 
origin in C and x and y axes along CD and CB, respectively. It is 0 otherwise. 
The double inequality at the beginning of this paragraph may be multiplied by 
x — xa and squared. A double inequality which is cubic in x is then obtained, 
from which the intervals in which the indicator function is 1 can be solved for 
and p conc bd',b'd(P) may be computed analytically. 

Pconc d'e'.de is a generic case and there is no simplification of functions 
(67) or (68). The expression for the measure density may be written as follows: 


Pconc D'E',DE\I) 



(80) 


where the integrand is a short hand notation for the function (69) showing the 
extremes of its indicator functions. The angle 7 is COD. 
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Figure 19: The five-pointed star 


11.2 The five-pointed star 

In Figure C3 the polar coordinates of the outer vertices of the five-pointed 
star are, starting from A, ( 1 , 7 r/10), ( 1 , 7 t/ 2 ), ( 1 , 97 t/ 10 ), (1, 137t/10) and 
(l,17vr/10). The polar coordinates of the inner vertices of the five-pointed 
star are, starting from B, (r, 37 t/ 10 ), (r, 77 t/ 10 ), (r, ll 7 r/ 10 ), (r, 37 t/ 2 ) and 
(r, 197 t/10), where r = *™( 7 n/io) • Due to the symmetry of the star, the set 
of chords which touch the side JA is congruent to the set of chords which 
touch any other side. If we multiplied the said set by 10, we would count each 
chord twice, because all chords (up to a set of measure 0) would touch two sides. 
Therefore the density of measure for the star is five times the the density of 
measure of the chords which touch the side JA. In its turn, 

PJA = PJA,AB + Pconc AJ,CD + PJA,DE + PJA,EF + PJA,FG■ (81) 

A rotation of 3(27t/5) transforms the pair of sides JA,DE into the pair of 
sides FG, JA. Therefore: 


where 


PJA — PJA,AB + Pconc AJ,CD + 2 PjA,DE + PJA,EF, 


(82) 


'\A-G\ 


Pconc AJ,Cd(I ) ~ dx ( 1 [GD,Gop(x)] (V-) - 1" 1 [GD,Gop(x)\ (V+) ' ' ' ) • (83) 

l\J-G | V 7 


The functions y± are given by formula (68) with 7 = 7 t /5 and the dots stand 
for 
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— xsiruTT/ j = ■ (arcsin (§ sin f) ± £), where the siqn of the last ~ matches 

the subindex of the argument of the corresponding indicator function. 

The above function is expensive to compute. The author has approximated 
the Riemann integral above by a Riemann sum with an interval of 0.00001. 
The computation of the graph that follows took 1,100 seconds in his personal 
computer, using the software Mathematica. Even after 1,100 one can still 
see some spikes due to numerical (not real) oscillations. The numerical na¬ 
ture of the spikes is revealed by the fact that they spread outside the bound 
[app-(n)(£),app + (n)(£)] (see section [2]). These spikes arise because for a given 
width of the Riemann sum the error is a rapidly varying function of £. 

R is numerically better to use the method given in section 0 When n = 
1000, the graphs of app-(n)(£) and app + (n )(£) for p conc aj,cd are indistin¬ 
guishable from each other. R takes 180 seconds to compute each of them, on 
the same computer as before. These bounds are exact in the sense that they 
are computed in terms of measure densities for convex polygons, which in turn 
we had computed in terms of elementary functions. 

We plot in Fig. 02| their average and the superposition with the graph 
Of pconc aj,cd{£ ) computed using the Riemann sum. Both ways of computing 
Pconc aj,cd are worked out in the file Concave.nb in ITW . 



Figure 20: Superposition of the graphs of p conc aj,cd{£ ) computed using (67) 
and (68) (which yields some small spikes) and using (39) 



Figure 21: Measure density of the length of chords for the five-pointed star 
In Figure HH the graph of the measure density (i. e., 5 Pja(£) ) is plotted. 
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12 Remark 

We notice that 



dx l \yi,y2\{y±(xJ)) 


x sm 7 


£\J £ 2 — x 2 sin 2 7 


sm 


^arcsin 


nr 

— sin 7 



[x 1 ,x 2 }n{x\y ± (x,t)&[y 1 ,y2]} 


x sin 'y / / x \ \ 

dx 7a o =?= sin arcsin 7 sin 7) ± 7) (84) 

£ \/£ 2 — a : 2 sm" 7 ' ' ' ' 


and that the primitives of 2 s i n (arcsin (f sin 7 ) ± 7 ) are 


sin 7 

~T 


x 


— ( x sin 7 T cot 7 


£ 2 


x 2 sin 2 7 


£ cos 7 

H-^— arctan 

2 sin 7 


\Jl 2 — x 2 sin 2 7 
x sin 7 


(85) 

The difficulty in computing the above integral then, reduces to the computa¬ 
tion of the integration limits. These limits are quite laborious to find, and have 
different analytical forms which depend on £ and on the geometry of the pair 
of segments. Still, it is possible in principle to obtain results in the manner of 
H. S. Harutyunyan, V. K. Ohanyan and A. Gaspary an f29i While this is 

less practical as a general purpose procedure to find distributions of chords, it 
may be worthwhile doing if one wishes to compute distributions of chords for 
a certain class of polygons repeatedly. 


Appendix: Opposite function 

We define the unit vectors 


v 


\ n 

\A — B 


and 


w 


C i) 
| C-D 


( 86 ) 


and let x and y be the distances to O along the straight lines OA and OC, 
respectively. 

We want to obtain the coordinates of O + yw when the coordinates of 
A, B , C, D and O + x v are known. Two ways to obtain them are the fol¬ 
lowing: Solve y from the equation 


(y w + O) — (x v + O) oc V — (x v + O) (87) 


or notice that the parallelograms which the vector O + x v — V spans with the 
vectors D — V and D — (O + y w) have equal areas, because their difference is 
parallel to O + xv — V. By either way one arrives at 
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Figure 22: The opposite function 


(O + x v - V) x (D - V) = (O + x v - V) x (D - O - y w), ( 88 ) 

from which y can be solved for: 


op(x) = y — \OD\ 


(O + x v - V) x (D - V) 
(O + x v — V) x w 


(89) 


When the segments AB and CD are parallel, 7 = BOD = 0, w = v and 
O lies at infinity. According to the previous scheme for 7 / 0 , the x and 
y coordinates would take infinite values at the segments. By convention, in 
this case, we take the origin of the x coordinates at B and the origin of the 
y coordinates at the closest point to B which lies on the straight line which 
contains the segment CD. As in the 7 ^ 0 case the chords which lie in the 
polygon under consideration are the ones which cross the triangle BVD. By 
the same argument as in the 7 / 0 case and with the said convention about 
the origin of coordinates x and y, 


op(x) 


y=(D - B).v 


{x v-(V- B)) x(D-V) 
(x v — (V — B)) x v 


(90) 
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